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Abstract 

The decay rates of P-wave heavy quarkonia to light hadrons are presented to leading order in v 2 and 
next-to-leading order in a s . They include contributions from both the color-singlet component and the 
color-octet component of quarkonia. Applying these results to charmonium and using measured decay 
rates for the Xci and Xc2 by E760, we determine the two nonperturbative decay matrix elements, and 
then predict the hadronic decay rates of XcO and h c , and the electromagnetic decay rates of XcO and Xc2- 
The obtained decay rates of XcO — * LH and XcO — > 77 are in agreement with the Crystal Ball result, and 
also with the new measurement by BES. However, the results for r(xcO — * LH) are dependent on the 
choice of renormalization scale. 
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The study of heavy quarkonium physics can provide very interesting tests of perturbative quantum 
chromodynamics (PQCD). Calculations of the rates for heavy quarkonium decay into light hadrons were 
among the early applications of PQCD. These early calculations are based on a naive factorization assumption 
that all long-distance nonperturbative effects can be factored into the nonrelativistic wavefunction of color 
singlet QQ or its derivative at the origin, and the perturbative part is related to the annihilation rates of 
color-singlet QQ which can be calculated using PQCD. In the nonrelativistic limit, this early factorization 
formalism was supported by explicit calculations for S-wave decays at next-to-leading order in a s But 
in the case of P-wave [g] quarkonium decays, infrared divergences appeared in the perturbative calculations 
of color-singlet QQ annihilation amplitudes. These are clear indications that the decay rates are sensitive to 
nonperturbative effects beyond those related to the wavefunction of color-singlet QQ pair or its derivative at 
the origin, and not all nonperturbative effects can be factored into the color-singlet component of quarkonium. 
Recently, Bodwin, Braaten and Lepage (BBL) have developed a rigorous factorization formalism || which 
is based on an effective field theory, nonrelativistic QCD (NRQCD). This factorization formalism provides a 
clean separation between short-distance effects and long-distance effects for the decay rates and production 
cross sections of heavy quarkonium. 

Nowadays there is a renewed interest in studying the decay of P-wave charmonium, not only due 
to the theoretical development mentioned above but also due to recent experimental results such as the 
total decay widths of Xcj and the observation of h c . BBL have applied the new factorization approach in a 
phenomenological analysis of P-wave charmonium decays Q . They give a leading order result with both the 
color-singlet and color-octet QQ components. The next-to-leading order correction to the decay of h c is given 
in pi, where both the color-singlet and color-octet contributions are included and the explicit cancellation 
of previously encountered infrared divergence is revealed. In the next-to-leading order color-singlet terms 
are considered in a phenomenological analysis of hadronic annihilation decays of XcJ- Recently the next-to- 
leading order color-octet corrections to hadronic \J decays have also been calculated [Q. In this paper we 
will perform a phenomenological study for the hadronic decays of four P-wave charmonium states by using 
the results that completely include the next-to-leading order QCD corrections. 

We start with the formulas for the P-wave quarkonium decay widths in the new factorization for- 
malism 

T( XJ -» LH) = 2Imf 1 ( 3 P J )H 1 + 2Imf s { 3 S 1 )H 8 + 0(v 2 T), (1) 

T(h - LH) = 2Imf 1 ( 1 Pi)H 1 + 2Imf 8 CS )H s + 0(v 2 T), (2) 

where Hi and Hg, are the matrix elements of color-singlet and color-octet operators respectively. The 
short-distance coefficients can be extracted by matching the imaginary part of the on-shell QQ pair forward 
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scattering amplitude calculated in full perturbative QCD with that calculated in NRQCD. We list the results 
to next-to-leading order in a s as the following 
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, where 

6 = i(llCU-2r*/), 

and Cf ~ ^iw^i @A — N c . The coefficients Imfi( 3 Pi) starts in order a 3 , hence they only are given to 
leading order; while all other coefficients start in order a 2 , whose next-to- leading order corrections are also 
given. We note that (||) and (|§|) are given in || and (Q) is given in Coefficients Imfx( 3 Pj) in (||)-(j|) have 
been calculated in ^| and listed in where quark and antiquark are taken off-shell and binding energy 
regularization scheme was used. Here we recalculate them by using dimensional regularization to control 
the infrared divergence, and there are some differences between our r esults (|)-(§) and that in @. This 
difference only comes from the diagram in Fig.l which represents the inclusive processes QQ( 3 Pj) — > qi<jig- 
In the d = 4 — 2e dimension space, contributions of the three particle cut diagram in Fig.l to the imaginary 
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part of QQ( 3 Pj) pair scattering amplitude are 
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Here 



(11) 



/(e) = (^rr(l + e) 

The results coming from the diagrams that represent the inclusive processes QQ( 3 Pj) — > g<? and QQ( 3 Pj) 
ggg are finite and have been given in p[. 



Fig.l Feynman diagram with three particle cut contributing to the divergence terms in the 
full theory calculation of QQ annihilation amplitudes 

While in the effective field theory NRQCD, the corresponding scattering amplitudes can be written 

as 

ImM{ 3 Pj) NRQCD = l^ll^A + (/ m/8( 3 5l))o i^(_ J_ + J_). (12 ) 

Comparing (|l2|) with (|9|)-(|ll|), it is obvious that the divergence terms are removed and finite coefficients 
Imfi( 3 Pj) (|J) — (||) can be obtained. It is important to point out that if one replaces In^ in the expressions 
in 0| by — 2e7^> then the divergent terms are the same as those in (^) — (pd|). The difference only occurs in 
their finite terms due to different regularization scheme being used. It is certainly true that the coefficients 
of 4-fermion operators must be infrared finite and independent of the choice of regularization procedures 
because all nonperturbative effects are factored into the matrix elements. Note that the coefficients can be 
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derived consistently only by taking the same regularization scheme in full QCD and in effective NRQCD. 
The advantage of using dimensional regularization is that the on-shell condition and gauge invariance are 
maintained manifestly and conventional treatment of NRQCD is under the on-shell condition, thus we 
can give an explicit cancellation for divergences appeared previously. The introduction of off-shell binding 
energy makes it difficult to do calculation in NRQCD, and the results are incomplete if simply absorbing the 
divergences associated with the logarithm of binding energy into the matrix elements of color-octet operators. 

Now we apply the factorization formula to charmonium systems. For the lowest radial excitation, 
the 3 Pj states are called Xcj and the 1 Pi state is called h c . The explicit form for their decay rates into light 
hadrons at leading order in v 2 are 

r( Xc0 -> LH) = C 00 a 2 (m c )(l + C i—)Hi + D a 2 s {m c )(l + D 1 —)Hs(m c )-, (13) 

TT TT 

Y{x* -> LH) = CicgHi + D a 2 (m c )(l + D x ^)H s {m c ), (14) 

TT 

T( X c2 -> LH) = C 20 a 2 s {m c ){\ + C 21 — + D a 2 s (m c ){l + D x — )H s (m c ), (15) 

TT TT 

T(h c -» LH) = CtfHi + D' a 2 s (m c )(l + D[^-)H s (m c ), (16) 

TT 

where "LH" on the left hand of (|l3|) — (|l^) represents all final states consisting of light hadrons, and the 
coefficients are 

4-7T 

Coo = C01 = 8.710; 
C\ = -0.370, C[ = -0.161; 

C 20 = ^, C 21 = -5.061; 
45 

D = tt, Dx = 4.110; 

In deriving these coefficients we have taken N c = 3, n/ = 3 and made a choice \x = m c for the scale in 

the MS scheme. The large size of some coefficients for the correction terms is apparent. These numbers 
obviously depend on the definition of the renormalizcd couplings a s . We will study the the renormalization 
scale dependence of the results later. In the following we use measured decay rates of the Xd and Xc2 to 
predict the inclusive decay rates of the Xco and h c , and the theoretical uncertainties will be estimated by 
considering relativistic corrections and high order perturbative QCD corrections. 
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Precision measurements of the total decay rates of the 3 Pi state Xai and 3 P 2 state %c2 have recently 
been carried out at Fermilab by the E760 collaboration. Their results with statistical and systematic errors 
are § 

r(x c2 ) = 2.00 ± O.lSMev, 
r(Xd) =0.88±0.14Mew. 

It is well known that the main decay modes of these P-wave charmonium states are the decay into light 
hadrons and the radiative transitions into J/ip or rj c . Other decay modes such as pionic transitions of 
the P states to the S states, of which the most important decay modes should be J/ip + 7T7t and r\ c + tttt, 
contribute much less to the total widths and therefore can be neglected ||. Previous experiments have 
measured the branching fractions for the radiative transitions of the Xci and Xci m to the J/ip, and they are 
B{Xd -> jJ/tp) = 0.273 ± 0.016, and B( Xc 2 -> jJ/tp) = 0.135 ± 0.011 §. We use the radiative branching 
fractions and the total decay rates to obtain the partial rates for light hadronic decays of the Xci and Xc2 

r( X ci -> LH) = 0.64 ± O.lOMev, (17) 

r(Xc2 -> LH) = 1.71 ± O.WMev. (18) 



Hi and H$ can be obtained directly by using (|14|) and (15) 



r (Xe2 lh) r( Xe i ^ lh) 

C 20 «2(l + C 21 ^)-C 1 a3 ' ^ y J 

Hs " ^fOTAf) ■ (20) 

Here we determine a B (m c ) by taking the coupling constant a s (mb) = 0.189±0.008 extracted from bottonium 
decays and evolving it down to the scale m c . The resulting value of the coupling constant is a s (m c ) = 
0.29 ± 0.02. Inserting (0) and © into © and ©, we obtain 

Hi = 18.4±5.2Mew, 

H R = 2.21 ± O.lbMev. 

The ratio of the two nonperturbative parameters is H%jH\ w 0.10, while it was determined to be 0.21 M if 



it was considered only to leading order in a s . Substituting Hi and H$ into ( |13|) and (16) we can easily get 
the decay widths of xm and h c into light hadrons 

r(x c o -> i-ff) = 12.4 ± 3.2Mev, (21) 



5 



T(h c -► LH) = 0.71 ± 0.07Mev. (22) 

Adding the radiative decay rate for Xco whose branching fraction has been measured to be (0.66 ± 0.18)%, 
we obtain the total decay rate r(xeo) = 12-5 ± 3.2Mev, which agrees with the earlier Crystal Ball value 
14 ± 5Mev ||, also with the new (preliminary) value of Ft ot (x c o) = 15.0^2'gMei> measured by BES, using 
3.5 x 10 6 -0' ("0(3686)) events in ip' — > 7XcO and Xco — * 7r + 7r~, K + K~ decay channels |l4|]. The rate for decay 
h c JVc has been estimated within a phenomenological framework to be about (340 — 380)kev |l3| ], and 
the total decay widths of h c is then T(h c ) « 1.07Mev, which is also consistent with the experimental result 

s 

For the electromagnetic decays to next to leading order in a s , we have 

r(XcO - 77) = 6^a 2 [l + (y - f )^]#i, (23) 

r(x c2 -77) = ye^ 2 (l-^)il 1 . (24) 
With the determined value for Hi, and e c = 2/3, a = 1/137, a s = 0.29, we predict 

r(xc0->77) = (3.72±l.ll)fcew 1 (25) 
which is in agreement with the observed value (4.0 ± 2.8)fcev by Crystal Ball Q, and also predict 

T(Xc2 -> 77) = (0.49 ± 0.W)kev, (26) 
which is larger than the E760 value [[ll] but smaller than the CLEO value |T^| . 

It is important to have a reasonable estimate for the theoretical uncertainties in our results. The two 
main sources of theoretical errors are relativistic corrections and higher-order perturbative corrections. Our 
formula (|l^)-(|l6|) are only valid to leading order in v 2 and high order relativistic corrections are not known 
at present. The error due to neglecting high order relativistic contributions could be of order v 2 rs 30%. On 
the other hand, we find that the one-loop coefficients in (jl3|)-(|l6|) are very large and strongly depend on the 
scale /x. It is well known that when working to all order in a s , the decay rates, being the physical observables, 
will not rely on the choice of /x. However we only do calculation to next-to- leading order in a s , therefore the 
analyses for the scale dependence of the results and the estimates for the higher-order effects are needed. 
The results are shown in Fig. 2 and Fig. 3 for decay rates r(x c o — > LH) and T(h c — > LH) respectively. For 

the running coupling constant a s with two-loops, three values = 200Gei>, 250Get», andSOOGev are 

used. The pictures show that our results are quite stable in the case of large /i, say, /i > 2m c . In the 
physically motivated range /i = m c to 2m c , the decay rates vary from 15Mev to 9Mev for T(xcO LH) and 
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Figure 1: Renormalization scale dependence of the decay width T(x c o — > LH) 

from 0.7Mev to 0.6Mev for T(h c — » LH) respectively, while the obtained two phenomenological parameters 
Hi = 22.0 - 19.0Meu, H s = 2.3 - 3.1Mev and the ratio H s /Ht = 0.11 - 0.16. It is interesting to note 
that, although there are large theoretical uncertainties due to the scale dependence and higher-order QCD 
corrections, our estimate for T(xcO ~ * LH) is enhanced greatly compared with the previous leading order 
result Ttot(Xco) = (4-8 ± 0.7)Mev jjj, which is smaller by a factor of three than the experimental value. 

Some comments on the relativistic corrections might be in order. It is obviously difficult to perform 
a complete analysis for the 0(v 2 ) corrections, because it must involve more higher dimensional four-fermion 
operators whose matrix elements are difficult to estimate at present. However, we might have some phe- 
nomenological analyses for the relativistic corrections in the color-singlet part. Just as the ratio T(xca ~^ 
2'y)/T(xc2 —> 27), which was discussed in ref.jy^, the color-singlet contribution to T(xcO — ¥ 2<?) /T(xc2 —> 2g) 
will receive relativistic corrections from two sources, i.e. the kinematic part and the dynamical part. In the 
language of the potential model, the color singlet matrix element Hi is proportional to R' p (0), the derivative 
of the wave function at the origin for the P-states. Due to a strong attractive spin-orbital force induced by 
one gluon exchange between quarks, which is also verified by the lattice calculations for the spin-dependent 
potentials between a heavy quark and an antiquark pq |, the Xco wave function in coordinate space will 
becomes narrower than the Xci wave function in which the spin-orbital force is repulsive, and therefore the 
derivative of the wave function at the origin becomes larger for Xco than that for Xc2- As a result, the 
dynamical relativistic effect will enhances Hi for Xco relative to Hi for Xc2, and this effect is found to be 
dominant over the kinamatic relativistic corrections 113] . This result might indicate that 0(v 2 ) corrections 
may further make T(xc0 ~ * LH) enhanced. As for the relativistic corrections in the color-octet part, more 
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Figure 2: Renormalization scale dependence of the decay width T(h c — > iff) 
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Figure 3: Renormalization scale dependence of H\ 
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considerations are apparently needed in the future work. 

In this paper we give the decay rates of four P-wave quarkonium states into light hadrons to leading 
order in v 2 and next-to-leading order in a s . They are expressed in terms of two nonperturbative parameters 
Hi and H$. Calculations in dimensional regularization scheme show that the infrared divergences, which 
appeared in the inclusive decay amplitudes for QQ — > qiqig and QQ — > ggg, can be cancelled explicitly by 
the contributions of olor-octet operators in NRQCD. The finite coefficients of Hi and Hs are given to next 
to leading order in a s . Using the derived theoretical results and the measured decay widths of Xci and Xci 
we estimate Hi, H% and the decay widths of Xco and h c . The determined values are very different from 
the previous values obtained by neglecting the next-to-leading order QCD corrections Q. In our results Hi 
is much larger than H%, and the decay width of Xco gets enhanced greatly due to 0{a s ) corrections. As a 
result, the predicted Xco hadronic decay width and electromagnetic decay width both could be in agreement 
with or close to the data. These significant differences indicate that QCD radiative corrections are very 
important in understanding the decays of P-wave quarkonium. However, our results are valid only to leading 
order in v 2 and next-to- leading order in a s . The large one-loop coefficients appearing in the expressions of 
the decay rates indicate that higher order QCD corrections may be important and that our results strongly 
depend on the choice of renormalization scale. More precise analyses must involve relativistic corrections 
and higher-order QCD corrections, which will include more matrix elements of higher dimensional operators. 

We would like to thank Professor E. Braaten for pointing out a numerical error in (^) and ( [TT] ) by 
comparing their recent result base on the threshold expansion method with our result by using the covariant 
projection method in dimensional regularization. It is turned out that the two methods in dimensional 
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regularization give identical results for the color-singlet sector of the P-wave decay widths, and are consistent 
with the previous calculation of Barbieri etal. using the binding energy as the infrared cutoff. 



References 

R. Barbieri, G. Curci, E. d'Emilio, and E. Remiddi, Nucl. Phys. B154,535(1979); K.Hagiwara, C.B.Kim, 
and T.Yoshino, Nucl. Phys. B177,461(1981). 

R. Barbieri, R. Gatto, and E. Remiddi, Phys. Lett. 61B, 465(1976); R. Barbieri, M. Caffo, and E. 
Remiddi, Nucl. Phys. B162, 220(1980); R.Barbieri etal., Nucl.Phys.B192,6l(198l). 

G. T. Bodwin, E. Braaten, G. P. Lepage, Phys. Rev. D51, 1125 (1995). 

G. T.Bodwin, E.Braaten, and G.P.Lepage, Phys. Rev. 46D, 1914 (1992). 



H.W.Huang, K.T.Chao, hep-ph/9601283, to appear in Phys. Rev. D 



M.L.Mangano and A.Petrelli, Phys. Lett. 352B, 445(1995) 
A.Pertrclli, CERN-TH/96-84 ( |hcp-ph/9603439| ) 



Particle Data Group, L. Montanet et al, Phys. #eu.D50(3-I),1171(1994). 

Y.P.Kuang, S.F.Tuan, and T.M.Yan, Phys. Rev. D37, (1988)1210. 

E760 Collaboration (T.A.Armstrong et al.), Phys. Rev. Lett. 69 (1992)2337. 

E760 Collaboration (T.A.Armstrong et al.), P/iys.i?e?;.Lett.68(1992)1468;70(1993)2988. 

CLEO Collaboration, J. Dominick et al., Phys. Rev. D50(1994)4265. 

K.T.Chao, Y.B.Ding, and D.H.Qin, Phys. Lett. B301, (1993)282. 

BES Collaboration, Y.F.Gu et al., presented at the Workshop on Beijing r — Charm Factory, Beijing, 
China, Feb. 1996. 

H.W.Huang, C.F.Qiao, and K.T.Chao, Phys. Rev D54, 2123(1996) 

A.Huntley, C.Michael, Nucl. Phys. B286, 211(1987); C.Michael, P.E.L.Rakow Nucl. Phys. B256, 
640(1985) 



10 



